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CN Abstract 

The behavior of a-quartz under hydrostatic and non-hydrostatic high-pressure conditions 
has been investigated in Molecular Dynamics (MD) simulations of silica in order to clarify 
the role of non-hydrostatic stresses in the amorphization process. It is shown that the amor- 
phization threshold is not modified if the stress along the c direction is lowered, so that the 
\q mean amorphization pressure can effectively be lowered under non-hydrostatic conditions. On 

the other hand, the application of a positive uniaxial stress along the c axis, results in the 
appearance of a new crystalline phase, where all silicon atoms are in five-fold coordination. 
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Previous experimental and theoretical studies have 
shown that a-quartz, the stable form of silica under 
standard conditions, undergoes a crystalline to amor- 
phous transition at high pressure and room tempera- 
ture [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. It has since been pro- 
posed that the application of pressure together with 
the inhibition of thermal activation led to a mechan- 
ical instability [6, 8] in the crystalline network driv- 
ing the crystalline to amorphous transition. It was 
proposed that non-hydrostatic stresses played an im- 
portant role in the amorphization process [4, 11] and 
it was shown in a recent experimental study on the 
high-pressure behavior of a-AlPO^ [12] that these 
stresses are of a fundamental importance in the crys- 
tal to amorphous phase transition of this material. 

In this letter, molecular dynamics (MD) simula- 
tions of the simple silica model recently developped 
in reference [13] are used to explore the effect of non- 
hydrostaticity on the crystal to amorphous transfor- 
mation. This interatomic potential [13], obtained 
by ab-initio calculations on H^SiO^ clusters, is cho- 
sen for its ability to reproduce a number of struc- 
tural properties of a-quartz, among which pressure 
induced amorphization at 22 GPa and 300 K [7]. 
The simulations are carried out in a modified (N,P,T) 
isobaric-isothermal canonical ensemble, in which the 
stresses cr xx ,cr yy ,cr zz can be varied independently. The 
simulation box is rectangular, with the z direction 
parallel to the c axis of a quartz and contains 1440 
atoms (480 SiO'i units). The external stress is mon- 
itored using Berendsen's scaling [14] independently 
in the three directions. The equations of motion are 
integrated with a time step of 1.34 fs. The densities 
of states reported in this paper are calculated from 
the (mass weighted) power spectrum of the velocity 
autocorrelation functions, computed in the (N,V,E) 
microcanonical ensemble. They are smoothed by a 
gaussian distribution of width 20 cm -1 , in order to 
have the same resolution as those obtained by inelas- 
tic neutron scattering experiements. 

A reference amorphous state was first created by 
high pressure amorphization under hydrostatic con- 
ditions. Amorphization occured at 22 GPa, as pre- 
viously reported [7, 8, 10]. That amorphous system 
was then compressed to 25 and 30 GPa. All three 
systems were allowed to age for 14 ps at constant P 
and T and for another 14 ps at constant V and T 
after which the vibrational densities of states (DOS) 
reported in figure 1 were computed. These systems 
were then decompressed to GPa during 130 ps and 
were aged during 14 ps at constant V and T once 
equilibrium was attained. The DOS were once again 
calculated and are also reported in figure 1. Co- 
ordination distribution before the transition (up to 
21 GPa) was identical to that of a-quartz, with all 



silicon atoms being four-coordinated with respect to 
oxygen. Table 1 lists the average Si-0 and O-Si coor- 
dination numbers above 22 GPa. The higher the pres- 
sure above the transition, the greater the proportion 
of six-coordinated silicons and three-coordinated oxy- 
gens. It can also be noticed that four bonded silicons 
and two bonded oxygens constitute the large majority 
of the decompressed amorphous material. This infor- 
mation is of prime interest because it is not possible 
to analyze SiO'i in- situ by X-ray absorption spec- 
troscopy due to the low energy of the silicon K-edge 
and high absorption of the diamonds in this energy 
domain. It is also clear that the higher the final pres- 
sure, the more the local environment upon decom- 
pression differs from the original material as far as co- 
ordination is concerned. The densities of state clearly 
reveal the growth of bands around 1000 cm -1 and 
between 600 and 700 cm -1 when pressure is raised 
after amorphization. The latter phenomenon has al- 
ready been observed by raman scattering on silica 
glass upon pressurization [15] and by infrared absorp- 
tion performed both on a-quartz and silica glass [3]. 
Upon decompression, the spectra reproduce qualita- 
tively the infrared measures previously reported by 
Williams et al. [3], at least in the high frequency 
range where the growth of the 900 cm -1 line can be 
observed. One should be very careful when analyz- 
ing the low frequency domain, since the finite size of 
the simulation box can affect the vibrational spec- 
trum. All together, these results confirm the remark- 
able ability of the model for describing silica under 
high pressure. 

Two different compression paths were then cho- 
sen for the application of non-hydrostatic stresses. 
First, the stress parallel to the c axis (cr zz ) was kept 
constant and below 22 GPa and the stresses in the 
perpendicular plane (cr xx and <T yy ) were raised up to 
amorphization. Then, a xx and a yy were maintained 
at a constant value below 22 GPa and cr zz was in- 
creased. 

In the first case, amorphization always occured 
for a xx and a yy between 22 and 23 GPa, hence the 
needed stress in the (a,b) plane for driving the phase 
transition does not seem to depend on the value of 
<r zz , at least in the range where the calculations were 
done. These runs were made with several values of 
a zz , namely 21, 20, 19, 18, 15 and 13 GPa. A close 
study of the radial distribution functions reveals that 
the very non-hydrostatic runs with low values of cr zz 
contained a small part of crystalline material. Si- 
O and O-Si coordination numbers were calculated 
in all six cases and are reported in table 1 for the 
most non-hydrostatic run. Upon decompression, all 
systems except the ones obtained under the most 
non-hydrostatic conditions remained amorphous, and 
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completely similar to those obtained under hydro- 
static conditions. The latter systems (a zz = lb and 
13 GPa) underwent a partially reversible transition 
in that sense that the recovered sample consisted 
of back-transformed quartz with some defects. The 
structure thus obtained is qualitatively similar to the 
one reported in figure 4. The Si-0 and O-Si coor- 
dination numbers as well as the vibrational density 
of state were calculated and compared with those of 
a-quartz (table 1 and figure 2). The first interpre- 
tation of these results for a-quartz, which is more 
compressible in the a and b directions than follow- 
ing the c axis, is that amorphization seems to hap- 
pen as the result of an instability occuring in this 
plane, or at least in a plane which has a small di- 
hedral angle with the latter [10, 11]. It also seems 
clear the the mean pressure of amorphization (i.e. 
the mean value of the three stresses) is lower than the 
pressure needed to hydrostatically amorphize the sys- 
tem. This is the first direct observation of this phe- 
nomenon, although most experimental studies seem 
to agree on this point. Amorphization is partially 
reversible upon decompression if the average non- 
hydrostatic pressure is low enough, and it could be 
that the crystalline domains at high pressure act as 
favorable nucleation sites upon decompression, driv- 
ing the amorphous-to-quartz transition. 

In the second series of runs, a xx and a yy were 
kept constant and a zz was raised until amorphization. 
This time, no direct quartz to amorphous phase tran- 
sition was observed. The system underwent a crystal- 
to-crystal transition towards a new phase, in which 
all silicon atoms are five-coordinated (see table 1). 
A slice parallel to the a-quartz c axis is represented 
in figure 3. The pressure conditions at which this 
structure appeared were a xx = a yy =20 GPa and 
<j zz =27 GPa. Pressure was isotropically raised from 
to 20 GPa, followed by an increase to 27 GPa along 
the c axis. Each compression was carried out over 
50000 steps accounting for a total compression time 
of 134 ps. 

In order to test the stability of this structure, its 
temperature was raised at constant pressure, with- 
out observable variation up to 1500 K. Upon further 
compression parallel to the c axis (cr zz =29 GPa), 
the sample became partially amorphous but in this 
case only five and six-coordinated silicon atoms were 
present in the disordered network (see table 1). 

The decompression phase was the exact opposite 
of the compression phase; the excess non-hydrostatic 
stress was first released and only then was the pres- 
sure lowered isotropically from 20 GPa. In the case 
of the new crystalline phase, a phase transition oc- 
cured at 2 GPa and the structure reverted back to the 
original quartz structure, whereas the partially amor- 



phous system transformed once again quasi-reversibly 
to quartz containing defect sites at 10 GPa (fig. 4). 
The corresponding DOS is reported in figure 2 and 
the coordination distribution in the table. 

This study clearly illustrates the importance of 
non-hydrostaticity in the solid state amorphization 
phenomenon. From a practical viewpoint, the most 
important conclusion is that the mean amorphization 
pressure can be considerably lowered in the presence 
of non-hydrostatic stresses. The amorphous solids 
that are produced under moderately non-hydrostatic 
conditions are very similar to those obtained in an hy- 
drostatic environment. The situation becomes quali- 
tatively different under strongly non-hydrostatic con- 
ditions (i.e. when the stress along the c axis differs by 
more than 30% of the other stresses). For small cr zz , 
amorphization is imperfect and the amorphous state 
reverts back to crystal upon decompression ("mem- 
ory glass"). It is interesting to note that the recov- 
ered material's vibrational density of state (fig. 2) 
presents broader spectral lines and weeker signal due 
to the apparition of defects upon recrystallization. 
These two observations have been reported experi- 
mentally in the case of a-AlPO^ berlinite [12] and 
anorthite [16] which both recrystallize upon decom- 
pression, and could therefore mean that these recov- 
ered crystalline phases contain some defects. For high 
<r zz , amorphization is pre-emted by a new crystal to 
crystal phase transition. The high pressure phase can 
amorphize, but again the resulting amorphous state 
is a "memory glass" that transforms back to quartz 
upon decompression. It has already been shown in a 
hydrostatic experiment [5] that phase transitions to a 
metastable crystalline form could constitute a path- 
way to solid state amorphization. This new phase, 
obtained under non-hydrostatic conditions could be 
another example of such crystal-to-crystal phase tran- 
sitions. 

These results are very suggestive of the great po- 
tential richness of solid state amorphization under 
non-hydrostatic conditions. It should be emphazised, 
however, that they were obtained using a very sim- 
ple model for silica, and very high compression rates. 
In particular, the existence of a new crystalline phase 
(which would be the first stable five-coordinated poly- 
morph of silica), and of reversible crystal to amor- 
phous transitions in silica, should be regarded as in- 
teresting possibilities, but are not established by such 
simulations. Carefuly designed non-hydrostatic anvil 
cell experiments will be necessary to investigate such 
effects. 
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Figure 1: The vibrational density of state upon compression and decompression in a hydrostatic environment. 
The pressures are reported on the figure . The spectra with the D superscript correspond to systems 
decompressed from the reported pressure. 
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Figure 2: The vibrational density of state of a system decompressed from a xx =a yy =22 GPa and cr zz = l3 GPa 
(thick solid line) and from a xx =a yy =20 GPa and a zz =29 GPa (thin solid line) compared with that of a- 
quartz at ambiant conditions (dotted line). The decompression in both cases is partly reversible. The 
sample back-transforms into crystalline quartz but contains some defects. Coordination of all three phases 
are almost the same (see table 1). This figure shows that the effect of defects present in the crystalline lattice 
is far from being negligible as far as the intensity and the width of the spectral lines are concerned. 



6 



Figure 3: The new high-pressure phase viewed down the c axis. The silicon atoms (large light grey spheres) 
positioned on the edges of the two similar equilateral triangles (4.5 angstroms of side length) are in the same 
plane and the distance between these triangles is apporximately 1 angstrom in the c direction. The oxygen 
atoms (small dark grey spheres) are disposed on the edges of two non-identical equilateral triangles. The 
atoms on the edges of these triangles constitue a sheet representing one third of the cell in the c direction. The 
two other sheets are superposed as to form a helix and the helical step is worth approximately 3 angstroms. 
Two sheets can be seen on this slice (indicated by the numbers near the silicon atoms). The equilateral 
triangle formed by the projection of the three helically superposed silicon atoms on the (a,b) plane has a 
side length of 1.1 angstroms. 
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Figure 4: A slice cut perpendicular to the c axis of the system decompressed from (20,20,29) GPa. The 
decompressed sample recrystallizes as quartz but contains some defects. The silicon atoms are represented 
as large light grey spheres whereas the oxygen atoms are the small dark grey spheres. This is a typical 
example of a "memory" glass. 
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Table 1: Si-0 and O-Si coordination distribution in hydrostatic and non-hydrostatic runs. The pressures 
followed by letter D correspond to a system at ambient pressure decompressed from the indicated pressure. 
In the first column, pressures represented with three figures in parentheses represent the <r xx , cr yy and a zz 
components of the stress, respectively. The last column lists the system's phase in the corresponding pressure 
conditions. 



P (GPa) 


Si IV 


Si v 


Si VI 


O n 


QUI 


Phase 





100% 


0% 


0% 


100% 


0% 


a-quartz 


22 


12.5% 


42% 


45.5% 


34.5% 


65.5% 


amorphous 


25 


8% 


38.5% 


53.5% 


30% 


70% 


amorphous 


30 


7% 


36.5% 


56.5% 


28% 


72% 


amorphous 


22 D 


75% 


23.5% 


1.5% 


87% 


13% 


amorphous 


25 D 


69% 


26% 


5% 


82% 


18% 


amorphous 


30 D 


56% 


36% 


8% 


73% 


27% 


amorphous 


(22,22,13) 


13% 


34% 


53% 


32% 


68% 


amorphous 


(22,22,13) D 


97% 


3% 


0% 


98% 


2% 


a-quartz + defects 


(20,20,27) 


0% 


100% 


0% 


50% 


50% 


new phase 


(20,20,27) D 


100% 


0% 


0% 


100% 


0% 


a-quartz 


(20,20,29) 


0% 


90% 


10% 


45% 


55% 


new phase + amorphous 


(20,20,29) D 


99% 


1% 


0% 


99% 


1% 


a-quartz + defects 
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